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New Results from the PIBETA Experiment 
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We report interim results of the PIBETA experiment analysis for the pion beta decay and pion radiative decay. The former is in excellent 
agreement with the SM predictions at the 1 % accuracy level. 



1 Experiment Goals and Motivation 

The PIBETA experiment['Tl at the Paul Scherrer Institute 
(PSI) is a comprehensive set of precision measurements of 
the rare decays of the pion and muon. The goals of the 
experiment's first phase are: 

(a) To improve the experimental precision of the pion 
beta decay rate, tt^ -n^e^v (known as TieT,, or nfi), 
from the present ~ 4 % to ~ 0.5 %. The improved 
experimental precision will begin to approach the 
theoretical accuracy in this process, and thus for the 
first time enable a meaningful extraction of the CKM 
parameter Yud from a non-baryonic process. 

(b) To measure the branching ratio (BR) of the radiative 
decay n — > evy {neiR, or RPD), enabling a precise 
determination of the pion form factor ratio FaIFv, 
and, hence, of the pion polarizability. Due to ex- 
panded phase space coverage of the new measure- 
ment, we also aim to resolve the longstanding open 
question of a nonzero tensor pion form factor. 

(c) A necessary part of the above program is an exten- 
sive measurement of the radiative muon decay rate, 
ju — > evvy, with broad phase space coverage. This 
new high-statistics data sample is conducive to a pre- 
cision search for non- (V-A) admixtures in the weak 
Lagrangian. 

(d) Both the nfi and the 7ii,2R decays are normalized to 
the TT — » ev (known as 11^2) decay rate. The first 
phase of the experiment has, thus, produced a large 
sample of Jiei decay events. The second phase of the 
PIBETA program will seek to improve the decay 
branching ratio precision from the current ~ 0.35 % 
to under 0.2 %, in order to provide a precise test of 
lepton universality, and thus of certain possible ex- 
tensions to the Standard Model (SM). 



Recent theoretical work[ |2| yj has demonstrated low the- 
oretical uncertainties in extracting Vud from the pion beta 
decay rate, i.e., a relative uncertainty of 5 x 10 or less, 
providing further impetus for continued efforts in improv- 
ing the experimental accuracy of this process. 



2 Experimental Method 

The ttEI beam line at PSI was tuned to deliver ~ 10* tt^/s 
with p„ ^ 113 MeV/c, that stop in a segmented plastic scin- 
tillator target (AT). The major detector systems are shown 
in a schematic drawing in Fig.[n Energetic charged decay 
products are tracked in a pair of thin concentric MWPC's 
and a thin 20-segment plastic scintillator barrel detector 
(PV). Both neutral and charged particles deposit most (or 
all) of their energy in a spherical electromagnetic shower 
calorimeter consisting of 240 elements made of pure Csl. 
The Csl radial thickness, 22 cm, corresponds to 12 Xq, and 
the calorimeter subtends a solid angle of about 80% of 
47r sr. 




Figure 1. A schematic cross section of the PIBETA detector sys- 
tem. Symbols denote: BC-thin upstream beam counter, AC1,2- 
active beam collimators, AD-active degrader, AT-active target, 
MWPCl,2-thin cylindrical wire chambers, PV-thin 20-segment 
plastic scintillator barrel. BC, ACl, AC2, AD and AT detectors 
are also made of plastic scintillator 

The basic principle of the measurement is to record all 
non-prompt large-energy (above the yU evv endpoint) 



2 



Workshop on the CKM Unitarity Triangle, IPPP Durham, April 2003 



electromagnetic shower pairs occurring in opposite detec- 
tor hemispheres (non-prompt two-arm events). In addition, 
we record a large prescaled sample of non-prompt single 
shower (one-arm) events. Using these minimum-bias sets, 
we extract nfi and 71^2 event sets, using the latter for branch- 
ing ratio normalization. In a stopped pion experiment these 
two channels have nearly the same detector acceptance, 
and have much of the systematics in common. 

A full complement of twelve fast analog triggers compris- 
ing all relevant logic combinations of one- or two-arm, 
low- or high calorimeter threshold , prompt and delayed 
(with respect to stop time), as well as a random and 
a three-arm trigger, were implemented in order to obtain 
maximally comprehensive and unbiased data samples. 
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Figure 4. Histogram of the measured y-y opening angle in pion 
beta decay events {n^ — » n'^e^v) for a subset of acquired data; 
solid curve: GEANT simulation. 
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Figure 2. rfi energy spectrum for a subset of the measured 
ifie^v decay data; solid curve: GEANT simulation. 
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Figure 3. Histogram of time differences between the beam pion 
stop and the nfi decay events (dots); curve: pion lifetime expo- 
nential curve. A software cut at 10 ns was applied. 



The high quality of the PIBETA data is clearly demon- 
strated in the histograms of the calorimeter energy and 
event timing (following the stop time), as well as of 
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Figure 5. Histogram of y-y time differences for the same set of nji 
data events (dots); curve: fit. Signal to background ratio exceeds 
250. 



the J -J opening angle and time difference for a subset of 
the recorded pion beta decay events, shown in Figs.|2j]5] 

In particular, the low level of accidental background is ev- 
ident in the y-y relative timing histogram in Fig. |5j the 
peak to background ratio exceeds 250. The histogram of 
recorded y-y opening angles for pion beta events, shown 
in Fig. 01 provides possibly the most sensitive test of the 
Monte Carlo simulation of the apparatus, and of the sys- 
tematics related to the geometry of the beam pion stopping 
distribution. The latter is the single largest contributor to 
the overall uncertainty in the acceptance, and, hence, in the 
branching ratio. 

3 First Results: Pion Beta Decay 

The first phase of measurements took place in 1999, 2000 
and 2001, resulting in some 60,000 recorded pion beta 
events. The figures of Section |2 are based on a data sub- 
set acquired in 1999 and 2000. Our current preliminary 
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working result for the pion beta decay branching ratio, ex- 
tracted from the above analysis, is 



BR ^ 1.044 + 0.007(stat.) + 0.009(syst.) x 10 



-8 



(1) 



Our result is to be compared with the previous most accu- 
rate measurement of McFarlane et al.[0: 



BR = 1.026 + 0.039 X 10"** , 

as well as with the SM Prediction (Particle Data Group, 
2002[|a): 

BR = 1.038 - 1.041 X 10-** (90%C.L.) 

(1.005 - 1.008 X 10"^ excl. rad. corr.) 

We see that even our working result strongly confirms the 
validity of the CVC hypothesis and SM radiative correc- 
tions [^,2j,3j|. Another interesting comparison is with the 
prediction based on the most accurate evaluation of the 
CKM matrix element Vud using the CVC hypothesis and 
the results of measurements of superallowed Fermi nuclear 
decays (Particle Data Group 2002[|5l): 

BR = 1.037 + 0.002 X 10"** . 

Thus, our current preliminary working result is in very 
good agreement with the predictions of the Standard Model 
and the CVC hypothesis. The quoted systematic uncertain- 
ties are being reduced as our analysis progresses. To put 
this result into broader perspective, we can compare the 
central value of Vud extracted from our data with that listed 
in PDG 2002[|5l: 



PDG 2002: 
PIBETA prelim: 



V,d = 0.9734(8), 
y„^ = 0.9771(56). 



Table [2 summarizes the main sources of uncertainties and 
gives their values both in the current analysis, and those 
that are expected to be reached in a full analysis of the en- 
tire dataset acquired to date. We have temporarily enlarged 
the systematic uncertainty quoted in Eq. ^pending a res- 
olution of the discrepancy found in the RPD channel and 
discussed in the following section. 

4 First Results: Radiative Pion Decay 

As was already pointed out, we have recorded a large data 
set of radiative decays: n:^ —> e^vy and e^vvy. To 

date we have analyzed both pion and muon radiative de- 
cays, though with more attention devoted to the former, as 
it is an important physics background to other decays un- 
der study. The pion radiative decay analysis has given us 
the most surprising result to date, and has commanded sig- 
nificant effort on our part to resolve the issue. 

The different event triggers used in our experiment are sen- 
sitive to three distinct regions in the RPD phase space: 



Table 1. Summary of the main sources of uncertainty in the ex- 
traction of the pion beta decay branching ratio. The column la- 
beled "current" corresponds to the present analysis based on a 
portion of the data taken. 

Summary of Dataset analyzed: 
uncertainties (%) partial* full 

external: 

pion lifetime 0.019 0.019 

BRin^ev) 0.33 -0.1^ 

BR{n° yy) 0.032 0.032 

internal: 

A(7T/3)/A(ev) 0.5 < 0.3 

At(y - e) 0.03 0.03 

E threshold < 0.1 < 0.1 

statistical: 0.7 ~ 0.4 

total: ~ 0.9 < 0.5 

* Subset of the 1999 and 2000 data. 
' Requires a new measurement. 



region A with e+ and y emitted into opposite hemi- 
spheres, each with energy exceeding that of the 
Michel edge (Em - 52MeV), recorded in the main 
two-arm trigger, 

region B with an energetic photon (Ey > Em), and 
Ee+ ^ 20 MeV, recorded in the one-arm trigger, and 

region C with an energetic positron (£'p+ > Em), and 
Ey ^ 20 MeV, also recorded in the one-arm trigger 



The RPD data are of a similar quality to our TTyS event set; 
due to limited space we omit the details here, and direct the 
interested reader to Ref. [ 1 1 instead. 

Together, the three regions overconstrain the Standard 
Model parameters describing the decay, and thus allow 
us to examine possible new information about the pion's 
hadronic structure, or non-(V-A) interactions. Appropri- 
ate analysis of these data is involved and nuanced, re- 
quiring a longer presentation than is possible here. We 
therefore only summarize the salient results of our work 
in progress on this pion decay channel. 

Our analysis indicates a measurable departure from SM 
predictions. Standard Model with the V-A electroweak 
sector requires only two form factors. Fa and Fy to de- 
scribe the so-called structure-dependent amplitude in RPD. 
The remainder of the decay amplitude is accounted for by 
QED in the inner-bremsstrahlung (IB) term. The pion vec- 
tor form factor is strongly constrained by the CVC hypoth- 
esis, while existing data on the radiative pion decay (PDG 
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2002[ 5 1) suggest that Fa - 0.5 Fy, yielding 

Fv = 0.0259 + 0.0005 , and Fa ^ 0.012 . 

Simultaneous as well as separate fits of our data in regions 
A, B and C confirm the above ratio of Fa/Fv - 0.5. How- 
ever, they show a statistically significant deficit in RPD 
yield in one region of phase space, for high Ey and lower 
Eg (mostly in region B), compared to predictions based on 
the above values of the pion form factors. 

A larger deficit in RPD yield, though less statistically sig- 
nificant than our result due to far fewer events, was first re- 
ported by the ISTRA collaboration[0|8l. This first obser- 
vation was interpreted by PoblaguevFlOl llOl as indicative of 
the presence of a tensor weak interaction in the pion, giving 

rise to a nonzero tensor pion form factor Fj 6 x 10"^. 

Subsequently, Peter Herczeg[ 11 1 found that the existing 
experimental evidence on beta decays could not rule out a 
small nonzero value of Ft of this order of magnitude. Ten- 
sor interaction of this magnitude could only be explained 
by the existence of leptoquarks. 
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Figure 6. Measured spectram of the kinematic variable A = 
(2Eg+/m^+) s'm^(9gy/2) in e^vy decay for the kinematic re- 

gion B, with limits noted in the figure. Solid curve: fit with the 
pion form factor Fv fixed by the CVC hypothesis. Ft = 0, and 
Fji free. Dashed curve: as above, but with Ft also released to 
vary freely, resulting in Ft = -0.0016 (3). Work in progress. 

We illustrate our working results in Fig. |6l which shows a 
projected one-dimensional distribution of A, a convenient 
kinematic variable based on Eg that ranges from to 1 re- 
gardless of Ey. It is clear that for lower values of A (and 
therefore of Ee), an SM fit with only Fv,Fa + overes- 
timates the experimental yield. Adding a nonzero tensor 
form factor of Ft 0.0016 produces statistically signif- 
icantly better agreement with the data. The fits are two- 
dimesional and encompass all three kinematic regions. A, 
B, and C. This work is in progress, and the reuslts are sub- 
ject to change — we are currently refining the analysis as 
well as the fit strategies. 



Taken at face value, this working result should not be in- 
terpreted as an indication of the existence of a tensor weak 
interaction, i.e., of leptoquarks. Instead, if it holds up in 
our final analysis, it would first suggest that the standard 
treatment of the RPD may not at this time correctly incor- 
porate all known SM physics. Radiative corrections seem 
to be a particularly good candidate for reexamination. 

5 Conclusions 

We have extracted an experimental branching ratio for the 
pion beta decay at the 1 % uncertainty level, and expect to 
reduce the uncertainty by another factor of about two in the 
near future. Our result agrees with the CVC hypothesis and 
radiative corrections for this process, and it opens the way 
for the first meaningful extraction of the CKM parameter 
Yud from a non-baryonic process. 

Our analysis of the n — > evy decay confirms that FaIFv - 
0.5, in agreement with the world average. However, events 
with a hard y and soft e^ are not well described by standard 
theory, requiring ''Ft + 0". A new theoretical look at this 
decay is needed. We can, though, rule out a large ''Ft" , as 
reported in analyses of the ISTRA data. 

The high statistics and broad coverage of our RPD data 
in principle guarantee extraction of pion weak form factor 
values with exceptionally low uncertainties. However, it 
appears that there may be significant theoretical uncertain- 
ties in the process of the form factor extraction. We hope 
that any remaining theoretical questions will be resolved in 
the near future, eventually enabling the field to make use 
of the full potential of the PIBETA data. 
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